Dielectric relaxation of DNA aqueous solutions 



S. Tomic,* S. Dolanski Babic,t and T. Vuletic 
Institut za fiziku, 10000 Zagreb, Croatia 

S. Krca and D. Ivankovic 

Institute Rudjer Boskovic, 10000 Zagreb, Croatia 

L. Griparic 

Department of Biological Chemistry, David Geffen School of Medicine, UCLA, Los Angeles, California 90095, USA 

R. Podgornik 

Department of Physics, University of Ljubljana and J. Stefan Institute, 1000 Ljubljana, Slovenia and 
Laboratory of Physical and Structural Biology NICHD, 
National Institutes of Health, Bethesda, Maryland 20892, USA 
(Dated: February 6, 2008) 

We report on a detailed characterization of complex dielectric response of Na-DNA aqueous 
solutions by means of low- frequency dielectric spectroscopy (40 Hz - 110 MHz). Results reveal 
two broad relaxation modes of strength 20 < Aelf < 100 and 5 < Aehf < 20, centered at 0.5 
kHz< !/LF <70 kHz and 0.1 MHz< ^-hf <15 MHz. The characteristic length scale of the LF process, 
50 < Llf < 750nm, scales with DNA concentration as Cp^'^^*""* and is independent of the ionic 
strength in the low added salt regime. Conversely, the measured length scale of the LF process 
does not vary with DNA concentration but depends on the ionic strength of the added salt as 
in the high added salt regime. On the other hand, the characteristic length scale of the HF 
process, 3 < Lhf < 50 nm, varyes with DNA concentration as c^^\ for intermediate and large 
DNA concentrations. At low DNA concentrations and in the low added salt limit the characteristic 
length scale of the HF process scales as cj^JJ,'^"^. We put these results in perspective regarding the 
integrity of the double stranded form of DNA at low salt conditions as well as regarding the role 
of different types of counterions in different regimes of dielectric dispersion. We argue that the free 
DNA counterions are primarily active in the HF relaxation, while the condensed counterions play a 
role only in the LF relaxation. We also suggest theoretical interpretations for all these length scales 
in the whole regime of DNA and salt concentrations and discuss their ramifications and limitations. 

PACS numbers: 82.39.Pj 87.15.He 77.22.Gm 



I. INTRODUCTION 



Semiflexible polyelectrolytes are ubiquitous in biolog- 
ical context ranging from charged biopolymers such as 
DNA or filamentous F actin, and then all the way to 
molecular aggregates such as bacterial fd viruses or the 
tobacco mosaic virus. They are an essential and fun- 
damental component of the cellular environment and 
make their mark in its every structural and functional 
aspect [1]. Their role is not confined solely to various 
(macro) molecular assemblies in the biological milieu but 
are equally prevalent in colloidal systems and soft mat- 
ter in general [2, 3]. It is their connectivity, stiiTness 
and strong electrostatic interactions, that allow polyelec- 
trolytes to show a wide range of complex behaviors, de- 
pending on their concentration, their overall length and 
the concentration and valency of the added salt ions and 
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intrinsic counterions [4]. 

In the simplest case of monovalent counterions, poly- 
electrolytes are usually stretched due to electrostatic re- 
pulsions and therefore statistically assume a rod-like con- 
figuration [5]. Polyvalent counterions on the other hand 
can turn electrostatic repulsion into attraction [6]. This 
s.c. correlation effect is one of the most important fea- 
tures of the polyvalent counterions and has fundamental 
repercussions for all charged soft matter. Correlation at- 
traction can strongly reduce the rigidity of charged poly- 
mers which then collapse into highly compact states [7] . 

Deoxyribonucleic acid (DNA) is in many respects a 
paradigm of a stiff, highly charged polymer. The struc- 
tural origin of the charge is due to negatively charged 
phosphate groups positioned along the DNA backbone 
[8] . The nominal charge density of double stranded DNA 
amounts to one negative elementary charge per 1.7 A or 
two elementary charges per base pair. When dealing with 
electrostatic interactions and their consequences in DNA 
one can usually ignore the internal chemical structure of 
DNA except as it transpires through the bare value of 
the persistence length (see below). 

Electrostatic interactions in stiff polyelectrolytes with 
monovalent counterions and added salt in aqueous solu- 
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tions arc standardly approached via the Poisson - Boltz- 
mann (PB) theory that combines electrostatics with sta- 
tistical mechanics on a simplified mean-field level [3] . On 
the PB level the outcome of the competition between en- 
tropy of mobile charges and their electrostatic interaction 
energy leads gross o modo to two types of effects. 

The first one is connected with electrostatic interaction 
of counterions and fixed charges on the polyelectrolyte. 
Positively charged counterions are attracted to the sur- 
face of the negatively charged DNA where they tend to 
accumulate in a layer of condensed counterions. This ac- 
cumulation can be understood within the framework of 
the Manning- Oosawa (MO) counterion condensation the- 
ory [3] as well as on the level of the solutions of the PB 
equation [9] . Within this theory the counterions accumu- 
late in the condensed layer around the cylindrical DNA 
surface only if the charge density parameter 77 > 1. Here 
T] = zls/b, where z is the valency of the counterion, b is 
the linear charge spacing and Ib is the Bjerrum length, 
defined as 



Ib = el/{4TveeokT) (1) 

Here eo is elementary charge, Eq is permeability of vac- 
uum, kT is thermal energy scale and e is the dielectric 
constant of the solvent. Bjerrum length is obviously de- 
fined as the separation between charges at which their 
electrostatic interaction energy equals their thermal en- 
ergy. In aqueous solutions = 7.1 A. The charge density 
parameter r] thus measures the relative strength of elec- 
trostatic interactions vs. thermal motion and is strongly 
dependent on the valency of the counterions. In the MO 
theory the counterions accumulate in the condensed layer 
exactly to such an extent that the effective charge den- 
sity parameter rj is reduced to 1 [10], i.e. the effective 
separation between charges is increased from b to Ib- 
The condensed counterions in the MO theory are still 
assumed to be perfectly mobile. Because of counterion 
condensation the effective charge of DNA is reduced by 
a factor r = 1 — l/{zr]). It is important to note here 
that the concept of counterion condensation is intrinsi- 
cally nonlinear and is a fundamental property of highly 
charged polymers. This includes DNA in both its double 
stranded {b = 1.7 A, rj = 4.2) as well as single stranded 
(6 = 4.3 A, ?7 = 1.7) forms. Counterion condensation 
for monovalent DNA salts is experimentally observed di- 
rectly by small- angle X-ray scattering measurements [11]. 
For double stranded DNA and monovalent counterions, 
counterion condensation occurs with r = 0.76. In the MO 
counterion condensation theory the condensation occurs 
only if the salt concentration is low enough to satisfy 
» a, where a is the polymer radius and is the 
Debyc screening length (see below) [10]. Furthermore the 
counterion condensation strictly occurs only in the limit 
of vanishingly small concentration of DNA [9] . For finite 
DNA concentrations a more complicated model of coun- 
terion condensation has to be invoked [9] that is based 



on the solution of the PB equation in the cell model ge- 
ometry. 

The second effect of competition between entropy of 
mobile charges and their electrostatic interaction energy 
is due to the interaction of mobile ions in solution be- 
tween themselves and their redistribution in the field of 
fixed polyelectrolyte charges. This redistribution leads 
to screening of electrostatic interactions between fixed 
charges. For small enough, fixed, charge density the 
screening is accurately described by the Debye-Hiickel 
equation, which is just a linearized form of the PB equa- 
tion [3] . On the Debye-Hiickel level the screening is quan- 
tified by a screening or Debye length defined as 

= SnlBn (2) 

where again Ib is the Bjerrum length and n is the den- 
sity of added salt. For monovalent salts the Debye length 

1/2 

is given numerically as 3.04/s , where Is is the 
ionic strength in M. Indeed for monovalent salts both ef- 
fects, ionic screening as well as counterion condensation, 
coexist. In general, however, one can not invoke screen- 
ing effects in polyvalent salt solutions since in that case 
the whole PB conceptual framework breaks down and 
correlation, not screening, effects [6] become the salient 
feature of the behavior of the system. The ionic screen- 
ing in monovalent salt solutions, augmented by the ef- 
fect of thermal DNA undulations, has been shown to 
quantitatively describe the measured osmotic pressure of 
DNA solutions in a fairly wide range of concentrations 
[12, 13, 14]. 

Apart from being a charged polymer, DNA is also 
molecularly rather stiff. The fiexibility of polymers is 
usually described via the persistence length Lp. Per- 
sistence length is nothing but the correlation length for 
orientational correlations along the molecular axis of the 
polymer [1]. For DNA the usually accepted value is about 
500 A [10]. Persistence length separates two regimes of 
behavior: rigid chain regime for contour lengths smaller 
than Lp and flexible chain regime for contour lengths 
much larger then Lp. The persistence length depends on 
long range electrostatic interactions along the polyelec- 
trolyte chain. The influence of electrostatic interactions 
on the persistence length was first considered by Odijk 
[15] and independently by Skolnick and Fixman [16]. Ac- 
cording to the Odijk-Skolnick-Fixman (OSF) theory, the 
total persistence length can be decomposed into a struc- 
tural (Lo) and electrostatic (ie) contribution as 

Lp = Lo + Le = Lo + lB/{2bKf (3) 

where b is again the separation between charges. As- 
suming the MO condensation, i.e. b = Ib, one gets 
Lp = Lo + 0.324/- 1 in A. As is clear from the OSF 
theory, counterion condensation reduces the electrostatic 

contribution to the persistence length due to an increase 
in the effective separation between charges from b to Ib- 
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The OSF result, though it can be nominaUy appHcd only 
at restrictive conditions, appears to work well when com- 
pared to experiments [17] as well as computer simulations 
[18]. On the other hand, in the regime of no added salt 
and thus weak electrostatic screening due to other chains 
and counterions, a semiflexible charged chain in a semidi- 
lute polyelectrolyte solution behaves like a random walk 
of correlation blobs (see below) with chain size R oc c~°'^^ 
[4]. Coming back to the electrostatic persistence length, 
it is worth mentioning that strongly depends on the 
valency of the counterions and in general on the details 
of the electrostatic interaction potential. For monovalent 
counterions Lg is usually positive indicating an effective 
repulsion between monomers. 

It is evident that starting assumptions ^ a of 
the MO theory (infinitely thin polymers or highly dilute 
polyelectrolyte solutions) and b of the OSF the- 

ory (monovalent salt with low enough screening) enforce 
their limited validity. It is not clear at all whether one can 
apply these theories in the limit where added salt concen- 
tration is smaller than the DNA concentration, viz. the 
concentration of intrinsic counterions. This fact raises 
further questions, even in the simplest case of DNA in 
monovalent salt solutions, of the general applicability of 
these models to predict the amount of counterion conden- 
sation, as well as to describe properties like fundamental 
length scales of DNA in solution. 

Due to the intrinsic length of native DNA, the concen- 
tration of polyelectrolyte solutions that we are dealing 
with in this contribution is always higher than the chain 
overlap concentration c* [4]. This means that we are ef- 
fectively always in the semidilute regime, where it has 
been known for a long time [4, 19], that a new length 
scale emerges describing the density correlations in the 
polyelectrolyte solution. This length scale is equal to the 
correlation length or the mesh size, describing the corre- 
lation volume in semidilute polyelectrolyte solutions, and 
is given by the de Gennes-Pfeuty-Dobrynin (dGPD) cor- 
relation length [4, 20, 21] that scales universally with the 
concentration of the polyelectrolyte, c, as 

e cx c-°-5 (4) 

It is interesting that this form of scaling, first derived 
for uncharged polymer solutions, is preserved also in the 
case of polyelectrolytes with some theoretically possible 
adjustments only in the case of added salt effects, and is 
expected to be proportional to the screening length due 
to both free DNA counterions and added salt ions [4] . For 
polyelectrolyte solutions the interpretation of this scaling 
result is that for volumes smaller than the polyelec- 
trolyte chain is stiffened by electrostatic interactions [22] , 
whereas for scales larger than ^, it behaves as a free flight 
chain. 

In an attempt to clarify several of the above issues, 
we have undertaken an investigation of dielectric relax- 
ation properties of DNA solutions that covers a broad 
range of DNA concentrations, as well as added salt 



concentrations. We have used low-frequency dielectric 
spectroscopy technique, widely established as a direct 
and non-destructive tool to probe charged entities and 
their structure in various bio-macromolecular systems 
[23, 24, 25]. Our aim was not only to verify the predic- 
tions of the various theories in the well defined conditions, 
but also to investigate the behavior of DNA at extreme 
conditions like very low and very high added salt limit. 
A brief report of this investigation has been published in 
Ref. [26]. 

Interpretation of our results depends on the nature of 
the conformation of DNA in low salt or even pure wa- 
ter solutions, more specifically on whether DNA is in a 
single stranded (ss-) or double stranded (ds-) conforma- 
tion. We have considered this question very carefully in 
what follows. Different preparation protocols for DNA 
solutions were adopted in order to study the issue of 
ds-DNA stability in pure water (see the Materials and 
Methods section) solutions. Denaturation of DNA was 
in particular investigated for pure water DNA solutions. 
These studies indicate that indeed for semidilute con- 
ditions the DNA double-helix was never denatured into 
two spatially distinguishable and well separated single 
strands. DNA solutions were additionally characterized 
by UV spectrophotometry, electrophoresis and atomic 
emmission spectroscopy. 

Our results demonstrate that DNA counterions, free 
as well as condensed in variable proportions, contribute 
to the oscillating polarization in the applied electric field 
and thus together determine the dielectric response of 
the DNA solution. The characteristic size of the relax- 
ation volume for the dielectric response of a semidilute 
DNA solution is given by one of the three fundamental 
length scales: the dGPD mesh size of the whole poly- 
electrolyte solution, the OSF salt-dependent persistence 
length of a single polyelectrolyte chain and the average 
size of the chain in the salt-free polyelectrolyte solution. 
The OSF prediction for the persistence length as a func- 
tion of added salt ionic strength is verified in the high 
added salt limit giving very good agreement. Our results 
indicate that by going from the high to the low salt limit 
and all the way down to the nominally pure water solu- 
tions, the characteristic length goes from the persistence 
length of the DNA in solution to a value that corresponds 
to the average size of the Gaussian chain composed of cor- 
relation blobs which scales as c^^'^^. Moreover, an exact 
condition is established that separates the high from the 
low added salt regime, which reads 21 s = 0.4ci„, where 
Cin is the concentration of DNA counterions, These two 
regimes differ in whether the added salt ions provide the 
screening or DNA acts as its own salt. Finally, our results 
confirm the theoretical prediction describing the concen- 
tration dependence of the mesh size of DNA solutions 
as ^ cx Cp^'^. In the limit of low DNA concentrations 
and low added salt, the semidilute solution correlation 
length deviates from the classical polyelectrolyte behav- 
ior and follows the scaling c^^'^''. Possible interpretations 
of this behavior are discussed and the appearance of lo- 
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cally fluctuating regions with exposed hydrophobic cores 
is suggested as the mechanism for this anomalous scahng. 

II. MATERIALS AND METHODS 

Salmon testes and calf thymus lyophilized Na-DNA 
threads were obtained from Sigma and Rockland, re- 
spectively. Low protein content was declared and ver- 
ified by our UV spectrophotometry measurements, for 
the former ^26oM280 = 1.65 — 1.70 and for the latter 
^26o/^280 = 1.87. Inductively coupled plasma - atomic 
emmission spectroscopy (ICP-AES) elemental analysis 
was performed on the Na-DNA threads dissolved in pure 
water [27]. The results have shown that the phospho- 
rus and sodium contents were the same, implying that 
only intrinsic sodium atoms are present. Taking into 
account that for double stranded DNA, monomers cor- 
respond to base pairs of molecular weight 660 g/mol, 
this result implies 7% by weight of Na"*" ions in DNA. 
This means that the concentration of intrinsic DNA 
counterions and the DNA concentration are related by 
Cin [mM] = cdna [mg/mL] x 3 /imol/mg. Electrophore- 
sis measurements were performed on DNA dissolved in 
pure water and in NaCl electrolyte with different ionic 
strengths. The obtained results have shown consistently 
the existence of polydisperse DNA fragments, most of 
them in the range 2-20 kbp. Since the scale of 3.4 A 
corresponds to one base pair, we estimate the range of 
contour lengths of DNA fragments to be 0.7 - 7 /xm. 

DNA solutions with different DNA concentrations and 
different added salt ionic strengths were prepared accord- 
ing to two protocols, which we describe below: 

I Pure water DNA solutions: Dry DNA threads 
were dissolved in pure water for 48 hours at 4^0 
so that the solutions within concentration range 
0.01 < Cdna < 15 mg/mL were obtained. The 
ionic strength of pure water Ig 0.01 mM was es- 
timated from the measured conductivity cr = 1.5 
/uS/cm of the pure water sample in the chamber 
for dielectric spectroscopy, using molar conductiv- 
ity of the highly diluted NaCl (126.5 Scm^/mol). 
An increased conductivity value of pure water, as 
compared to the declared one [27], is due to manip- 
ulation in the laboratory environment. 

II DNA solutions with added salt: 

111 NaCl was added to DNA water solution with a 
chosen cdna (prepared according to I), so that 
the added salt ionic strength was achieved in 
the range 0.01 mM< < 4 mM. 

112 DNA solutions with the same cdna as in III 
and with the added salt ionic strength in the 
range 1 mM< /« < 4 mM were prepared start- 
ing from stock DNA solutions in which DNA 
was dissolved in 10 mM NaCl for 48 hours at 
4°C. One of the stock solutions was dialyzed 



against lOmM NaCl during 24 hours at 4°C 
(112. 2), while the other was not (112. 1). 

113 DNA solutions with concentrations in the 
range 0.1 < cdna < 1-25 mg/mL and with 
the added salt ionic strength /« = 1 mM were 
prepared starting from a stock DNA solution 
in which DNA was dissolved in 10 mM NaCl 
for 48 hours at 4°C. Stock solution was dia- 
lyzed against 1 mM NaCl during 24 hours at 

rc. 

The stock solutions were stored at 4°C and were di- 
luted just before dielectric spectroscopy measurements, 
which were completed in a few days. UV spectrophotom- 
etry and electrophoresis measurements were performed 
within next two weeks. For longer periods the stock solu- 
tions were stored at -80°C. pH of all solutions were found 
to be around 7. Similar pH values were measured for 
DNA solutions prepared from 10 mM Tris-EDTA buffer 
(TE), which adjusts pH to 7.6 (not discussed in this ar- 
ticle). 

The results of ICP-AES elemental analysis indicated 
that solutions of Na-DNA dissolved in pure water did not 
contain any additional ions except declared sodium ones. 
This was further confirmed by dielectric spectroscopy 
measurements, which gave identical results for DNA solu- 
tions with same DNA concentration and added salt ionic 
strength prepared according to the protocol III and 112. 

Spectrophotometry measurements at 260 nm were per- 
formed to verify nominal DNA concentrations. The con- 
centration was determined assuming an extinction coef- 
ficient at 260 nm A260 = 20 for 1 mg/mL, meaning that 
the measured absorption A2eo = 1 corresponds to 0.05 
mg/mL [28, 29]. Throughout this article we will refer 
to nominal concentrations, which we have found to be 
in a good agreement with the measured ones. In par- 
ticular, the measured concentrations for DNA solutions 
112 and 113, performed on 10 aliquots, were consistently 
smaller by about 20% than the nominal ones. We in- 
terpret this difference to be due to water content, not 
taken into account by the spectrophotometry approach. 
Indeed, lyophilized DNA if kept at 110°C for 30 min loses 
about 20% in weight. 

Dielectric spectroscopy measurements [25] were per- 
formed at room temperature (25°C) using a set-up which 
consists of a home-made parallel platinum plate capaci- 
tive chamber and temperature control unit, in conjunc- 
tion with the Agilent 4294A precision impedance ana- 
lyzer operating in = 40 Hz - 110 MHz frequency range. 
The capacitive chamber enables reliable complex admit- 
tance measurements with reproducibility of 1.5% of sam- 
ples in solution with small volume of 100 fiL and with 
conductivities in the range of 1.5 - 2000 ^S/cm. Low ac 
amplitudes of 50 mV were employed in order to probe 
the DNA response in the linear regime, once we verified 
that for ac signal levels in the range between 20 mV and 
500 mV the result was essentially the same. Lowest ac 
amplitudes were not used in order to avoid extensive aver- 
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aging. Admittance was sampled at 201 frequencies at 27 
points per frequency decade. At each frequency, admit- 
tance was sampled 10 times and averaged. In addition, 
three consecutive frequency sweeps were taken in order 
to average out the temperature variations. Total time for 
described measurement amounts to 60 sec. Measurement 
functions are the real part of the admittance Gexp(i^) and 
the capacitance Cexpi^^), where ui = 2'ku. 

In addition to DNA samples, the reference sam- 
ples were also measured in order to minimize stray 
impedances, including the free ion contribution and elec- 
trode polarization effects, and extract the response due 
to DNA only [25, 30]. Reference samples were cho- 
sen as NaCl solutions of different molarities, adjusted 
to have the real part of admittances at 100 kHz the 
same as DNA solutions. In comparison to these solu- 
tions, independent adjustment of capacitance at 1 kHz 
demanded the reference NaCl concentrations that differ 
at most by 20%. Since the difference in capacitance be- 
tween reference NaCl solutions of similar molarities is 
approximately constant at frequencies above the influ- 
ence of the electrode polarization, adjusting solely the 
real part of admittance effectively also adjusts the ca- 
pacitance up to an additive term. As a result, the sub- 
traction of the reference response successfully eliminated 
influence of spurious effects down to a low frequency 
limit in the range 0.5 - 30 kHz, depending on the so- 
lution molarity, and up to a high frequency limit of 30 
MHz. Generally speaking, electrode polarization effects 
are larger for higher DNA and added salt concentrations, 
so that the low frequency bound is shifted to higher fre- 
quencies. Assuming that the conductivity contribution 
of each entity in the solution is additive [2], the DNA 
response is given by G{uj) = Gexp{^) — Gref{^) and 

C{uj) = Cexp{^) - Cref ioj), where Gre/(w), Cref{oj) is 

the reference samples response. Finally, the real and 
imaginary parts of dieletric function are extracted using 
relations 



phenomenological Havriliak-Negami (HN) function [31] 

1 



e{u!) — shf = Ac 



I C{lo) 
S so 
I G{uj) 
S ioso 



(5) 

(6) 



l/S = 0.1042 ± 0.0008 cm-i is the chamber con- 
stant, where S = 0.98 cm^ is the effective electrode 
cross-section corresponding to the sample of 100 /iL and 
I = 0.1021 ± 0.0001 cm is the distance between the elec- 
trodes. The chamber constant was determined by mea- 
suring the real part of admittance at 100 kHz of a 0.01 M 
and 0.0005 M KCl standard solutions (Mettler- Toledo) . 
It was also corroborated by the difference in capacitance 
of the chamber when empty and with 100 /uL of the pure 
water. In the latter case, the dielectric constant of pure 
water = 78.65 was used as a standard. 

Detailed analysis of the DNA response was made in 
terms of the complex dielectric function e{uj) given by 
a generalization of the Debye expression known as the 



{1 + {t^'Tay-'-f 



(7) 



where Ae = eq — ehf is the strength of the relaxation 
process, Eq is the static dielectric constant {u <C 1/to) 
and £hf is the high frequency dielectric constant {lo ^ 
1/to). To is the mean relaxation time, while 1 — a and /? 
are the shape parameters which describe the symmetric 
broadening of the relaxation time distribution function 
and skewness, respectively. We worked with f3 = 1 since 
this simplified HN formulation (also known as the Cole- 
Cole function) has been widely and successfully used to 
describe relaxation processes in disordered systems. 

Measured data were analyzed by using the least 
squares method in the complex plane [32, 33]. Such 
an approach, which takes into consideration both the 
real and the imaginary part of the dielectric function at 
the same time, strongly improves the resolution if com- 
pared with the method in which the real and imaginary 
parts are treated separately. The complex plane method 
proved itself to be a powerful tool to resolve reliably two 
close modes in frequency even if the strength of one or 
both modes docs not exceed Ae = 2 — 3, providc^d that 
the ratio of their dielectric strengths is smaller than the 
ratio of their positions in frequency. 



III. RESULTS 

We now present results of the dielectric response study 
performed on DNA solutions prepared according to pro- 
tocols described in Section II. The obtained dielectric re- 
laxation results for salmon testes and calf thymus DNA 
samples were essentially the same. For the sake of sim- 
plicity, in the remainder of this article we will speak of 
DNA samples only. 

Fig. 1 shows the frequency dependent real and imag- 
inary part of the dielectric function for selected DNA 
concentrations of DNA solutions. The results for pure 
water DNA solutions (protocol I, concentrations al=2.5 
mg/mL, a2=0.4 mg/mL, a3=0.1 mg/mL and a4=0.0125 
mg/mL) are shown in panel a) and b), while the results 
for DNA solutions with added salt of ionic strength Is = 
1 mM (protocol 113, concentrations bl=0.83 mg/mL, 
b2=0.5 mg/mL, b3=0.31 mg/mL and b4=0.125 mg/mL) 
arc shown in panel c) and d) . The observed dielectric re- 
sponse is complex and the data were only successfully 
fitted to a formula representing the sum of two HN func- 
tions. The full lines in Fig. 1 correspond to these fits, 
while the dashed lines represent single HN forms. The 
main features of this response, for pure water DNA so- 
lutions, as well as for DNA solutions with added salt, 
arc two broad modes, whose amplitude and position in 
frequency depend on the DNA concentration. The pa- 
rameter 1 — a, which describes the symmetrical broad- 
ening of the relaxation time distribution function, is 
concentration independent and similar for both modes 
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FIG. 1: Double logarithmic plot of the frequency dependence 
of the real (e') and imaginary (e") part of the dielectric func- 
tion at T = 25°C of (a, b) pure water DNA solutions (protocol 
I) and (c, d) DNA water solutions with added salt Is = 1 mM 
(protocol 113) for representative al-a4 (2.5, 0.4, 0.1, 0.0125 
mg/mL) and bl-b4 (0.83, 0.5, 0.31, 0.125 mg/mL) DNA con- 
centrations. The full lines are fits to the sum of the two HN 
forms; the dashed lines represent a single HN form. 



1 ^ a 0.8. The mode centered at higher frequencies 
(0.1 MHz< fHF <15 MHz) is characterized by smaller 
dielectric strength (5 < Achf < 20) than the mode 
(20 < Aelf < 100) centered at lower frequencies (0.5 
kHz< j^LF <70 kHz). In the remainder of this article, we 
will refer to these modes as the high-frequency (HF) and 
low-frequency (LF) mode, respectively. 

In what follows, we discuss possible assignments for 
these relaxation modes inside the framework of exist- 
ing theoretical approaches for polarization response of 
charged biopolymers in solution [25, 35, 36, 37, 38]. An 
applied ac field generates an oscillating fiow of net charge 
associated with DNA counterions [39] and induces polar- 
ization. Since the counterion displacement is controlled 
by diffusion, the dielectric response is basically charac- 
terized by the mean relaxation time tq cx L'^ /Din, where 
L is the associated length scale, and Din is the diffusion 
constant of counterions which is sufficiently well approx- 
imated by the diffusion constant of bulk ions [25, 40]. 
Since wc deal with Na-DNA solutions, we take the diffu- 
sion constant of Na+ ions Din = 1-33 • 10~^ m^/s. Note 
that the equation tq oc L"^ /Din is a scaling relationship 
and the proportionality constant is of order one (see Dis- 
cussion, part B). 

Several length scales are theoretically expected to be 
associated with dielectric relaxations of polyelectrolytes 
in solution: the contour length, the Debye screening 
length, the polymer chain statistical segment length and 
the polymer solution mesh size. 

The mean relaxation time for pure water DNA solu- 
tions, as well as for DNA solutions with added salt, is 
found in the range 10~^ — 1.5 • 10~^ s for the HF mode 
and 2 • 10"^ - 3 • 10"'' s for the LF mode. The correspond- 
ing characteristic length for the HF mode spans the range 
from 3 nm to 50 nm, while the characteristic length for 
the LF mode varies between 50 nm and 750 nm. Both 
of them are thus not within the range of the contour 
length distribution in our samples. On the other hand, 
the Debye screening length, the polymer solution mesh 
size and the polymer statistical segment length appear 
as plausible candidates. Moreover, it is noteworthy that 
specifically the values of the characteristic length for the 
LF mode are close to the values expected for the DNA 
persistence length. 



A. HF mode 

First, we address the HF mode. In the case of pure 
water DNA solutions the characteristic length Lhf shows 
the scaling Lhf oc c^^'^ with respect to the DNA con- 
centration, all the way down to a crossover concentration 
Ceo ~ 0.6 mg/mL. At that point the scaling form is then 
changed to Lhf « c^na'' (Fig- 2a). The observed behav- 
ior at high DNA concetrations conforms exactly to the 
de Gennes-Pfeuty-Dobrynin (dGPD) [4, 20, 21] scaling 
form valid for salt-free polyelectrolyte solutions [4, 25]. 
At low DNA concetrations, in the regime below Cco, -^^hf 
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FIG. 2: a) Characteristic length of the HF mode (Lhf) for 
pure water DNA solutions (protocol I, open squares) and for 
DNA solutions with added salt Is — 1 mM (protocol 113, full 
squares) as a function of DNA concentration (cdna). The 
full line is a fit to the power law Lhf occ^^'^^ and occ^^'^ for 
Cdna smaller and larger than Cco ~ 0.6 mg/mL, respectively, 
b) Characteristic length of the HF mode (Lhf) for DNA so- 
lutions with varying added salt (Is) for three representative 
DNA concentrations: cdna ~ 0.05 mg/niL (diamonds, pro- 
tocol III), Cdna = 0.5 mg/mL (full triangles, protocol III; 
open triangles, protocol 112. 1; open inverse triangles, proto- 
col 112. 2) and cdna = 0.8 mg/mL (circles, protocol 112. 1). 
The full line denotes Debye screening length for the in- 
vestigated range of added salt ionic strength 7s. 



displays an unusual behavior generally not observed in 

scmidilutc solutions. It appears as though at low DNA 
concentrations local conformational fluctuations partially 
expose the hydrophobic core of DNA so that the corre- 
lation length scales as Cj^^'^'^, as is the case of charged 
chains with partially exposed hydrophobic cores [4]. In 
dilute solutions, but only in dilute solutions, this scaling 
form would be typical for the average separation between 
chains [4]. The observed crossover might be thought to 
reflect the border between dilute and semidilute solutions 
[25, 41] corresponding to the crossover concentration c* 
[20] . For the shortest fragments of 2 kbp in our DNA so- 
lutions, we get c* of the order of 0.006 mg/mL, while the 
lowest concentration of DNA solutions is 0.01 mg/mL. In 



this manner, the interpretation of Cco ~ 0.6 mg/mL as 
the dilute-semidilute crossover concentration c* is ruled 
out. 

With added 1 mM salt, the dGPD behavior of Lhf 
remains unchanged, thus Lhf oc Cj^j^'^, as long as the 
concentration of intrinsic counterions Ci„ (proportional 
to Cdna) is larger than the concentration of added salt 
ions 21 s (Fig. 2a). When the concentration of intrinsic 
counterions becomes smaller than the added salt concen- 
trations, the Lhf apparently shows a leveling off , with 
a limiting value close to the Debye length appropriate 
for this salt concentration. One should be cautious here 
since the data become much less reliable exactly at low 
DNA concentrations (see error bars in Fig. 2a). Three 
sets of additional data (Fig. 2b) for three representative 
DNA concentrations with varying added salt also seem to 
reveal that Lhf does not vary with Ig in most of the mea- 
sured range of added salt, while the corresponding Debye 
screening length in the same range of added salt values 
decreases substantially. Unfortunately the accuracy of 
the data again becomes much less reliable due to the 
progressive merging of the HF and LF modes when one 
approaches the regime, where the characteristic length 
scale becomes apparently larger than the nominal Debye 
length at that salt concentration. 

Next we consider the behavior of dielectric strength 
defined as Aehf ~ /hf • Cm ■ cchf, where /hf is the 
fraction of counterions participating in the HF process, 
Cin [mM] = Cdna [mg/mL] x 3 /imol/mg (as explained 
in Section II) and ckhf is the corresponding polarizabil- 
ity. The polarizability ckhf is given by the scaling form 
aHF oc e^-L^p-/(eo^r) oc Ib-s-L'^^, [25, 41]. Therefore, 
the fraction of counterions /hf participating in the HF 
process is proportional to A£hf/(cdna • ^hf)- I^iS- ^ 
we show the dependence of /hf on DNA concentration 
(panel a)) and on the ionic strength of added salt ions 
(panel b)). 

The /hf data for pure water DNA solution displayed in 
panel a) indicate that the fraction of counterions partici- 
pating in this relaxation process does not depend on the 
concentration of DNA. Since the HF relaxation happens 
at the length scale ^ which describes the density corre- 
lations between DNA chains, this result indicates that it 
is the free counterions as opposed to condensed counteri- 
ons, that can hop from chain to chain in the volume 
that are the relaxation entities participating in the HF 
process. It is noteworthy that a similar interpretation 
was previously proposed by Ito et al. [41] for the relax- 
ation in synthetic polyelectrolytes observed in the same 
frequency range [25]. 

The data displayed in Fig. 3 b) suggest that the frac- 
tion of intrinsic counterions /hf active in the HF mode 
remains constant when salt is added to the DNA solution 
as long as cdna is substantially larger than L,. However, 
as soon as the concentration of added salt ions prevails 
over the concentration of intrinsic counterions, /hf starts 
to decrease. This is also discernible in the behavior of 
the /hf data for Is = 1 mM added salt solution shown 
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FIG. 3: a) Normalized dielectric strength of the HF mode 
AeHF/(cDNA • inp) as a function of DNA concentration 
(cdna) for pure water DNA solutions (protocol I, open 
squares) and for DNA solutions with added salt Is = 1 mM 
(protocol 113, full squares). The full line is a guide for 
the eye. b) Normalized dielectric strength of the HF mode 
AeHF/(cDNA • ier) vs. ionic strength of the added salt 
(Is) for three representative DNA concentrations: cdna = 
0.05 mg/mL (diamonds, protocol III), cdna = 0.5 mg/mL 
(full triangles, protocol III; open triangles, protocol 112. 1; 
open inverse triangles, protocol 112. 2) and cdna = 0.8 mg/mL 
(circles, protocol 112. 1). 



in panel a). A plausible suggestion would be that the 
salt renormalization of /hf is a consequence of screening 
due to added salt ions that seem to diminish the effective 
number of counterions that can participate in the chain 
- chain hopping process. 



B. LF mode 

Second, we address the LF mode. For pure water DNA 

solutions (protocol I), the characteristic length Llf in- 
creases with decreasing DNA concentration in almost 
three decades wide concentration range (open inverse 
triangles in Fig. 4a)) following the power law Llf oc 
Cdna The exponent -0.29 ± 0.04 suggests that 



in this regime _Llf is proportional to the average size of 
the polyelectrolyte chain that behaves as a random walk 
of correlation blobs and scales as c^^'^^ [4]. For DNA 
solutions with added salt /g = 1 mM (protocol 113, full 
inverse triangles in Fig. 4a)), Llf coincides with the one 
found for pure water solutions with high DNA concentra- 
tions. As soon as the concentration of intrinsic counteri- 
ons Cin (proportional to Cdna) becomes smaller than the 
concentration of bulk ions from added salt 2Is , Llf starts 
to deviate from the Llf oc Cj^na^ behavior and decreases 
to attain value of about 500 A at which it saturates. 

The dependence of Llf on the added salt ionic strength 
Is is shown in Fig. 4b) for three DNA concentrations. 
The observed data can be nicely fit to the OSF behavior 
[15, 16] of the form: Lp = Lo + a- 1''^. We get Lq = 470 
A for the structural persistence length and a = 0.09 
Am. While the value of Lq close to 500 A is in accordance 
with standard expectations for DNA [10], the value of the 
coefficient a is somewhat smaller than expected by the 
OSF theory a = 0.324 AM. It is noteworthy that the OSF 
model applies as long as the ionic strength of added salt is 
larger than the concentration of the intrinsic counterions. 
The data for cdna = 0.05, 0.5 and 0.8 mg/mL deviate 
from the OSF behavior for 1^ < 0.03 mM, < 0.3 mM 
and 7s < 0.5 mM, respectively. The value of Llf in this 
low salt limit attains the same value as in pure water 
DNA solutions (see Fig. 4a). 

In Fig. 5a) we show dependence of /lf oc Aelf / (cdna • 
on the DNA concentration assuming again that 
the polarizability Olf varies as the characteristic length 
squared. The data show that the fraction of counterions 
/lf active in the LF mode is roughly speaking indepen- 
dent of the DNA concentration. This is valid for the pure 
water DNA solution, as well as in the case of Is = 1 mM 
added salt solution. In our view this result together with 
the fact that the LF relaxation happens at the length 
scale of the average size of the polyelectrolyte chain sug- 
gests that the LF relaxation engages mostly condensed 
counterions along and in close vicinity of the chain. 

However, the data displayed in panel b), as well as 
the data for Ig = 1 mM in panel a) suggest that the 
fraction of counterions participating in the LF process 
/lf becomes larger in the case of added salt solutions, 
compared to the pure water case, if the concentration of 
added salt ions becomes larger than cdna- This means 
that at least some of the free counterions join the relax- 
ation of the condensed counterions along the segments of 
the same chain. It is thus impossible to completely sepa- 
rate condensed counterions from free counterions in their 
contribution to the LF relaxation mode. This conclusion 
bears crucially on the assumption that the scale of polar- 
izability is given by Llf(7s). It is noteworthy that for the 
HF and LF modes the addition of salt changes the effec- 
tive number of participating counterions in the opposite 
way. This might be attributed to the increased screen- 
ing for the interchain HF relaxation, and to the intrinsic 
counterion atmospheres squeezed closer to the chains due 
to reduced Debye length for the LF relaxation along the 
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FIG. 4: a) Characteristic length of the LF mode (Llf) for 
pure water DNA solutions (protocol I, open inverse triangles) 
and for DNA solutions with added salt Is = 1 niM (pro- 
tocol 113, full inverse triangles) as a function of DNA con- 
centration (cdna)- The full line is a fit to the power law 
I/LF oc Cq^'^^*'' ° . b) Characteristic length of the LF mode 
(Llf) for DNA solutions with varying added salt (7s) for three 
representative DNA concentrations: cdna = 0.05 mg/mL (di- 
amonds, protocol III), Cdna ~ 0.5 mg/mL (full triangles, pro- 
tocol III; open triangles, protocol 112. 1; open inverse trian- 
gles, protocol 112. 2) and cdna = 0.8 mg/mL (circles, protocol 
112. 1). The full line is a fit to the expression Lp = Lo + a- 1^^ 
with Lo = 470 A and o = 0.09 A. 



FIG. 5: a) Normalized dielectric strength of the LF mode 
Aelf /(cdna -I/lf) b£ a function of DNA concentration (cdna) 
for pure water DNA solutions (protocol I, open inverse tri- 
angles) and for DNA solutions with added salt 7s = 1 mM 
(protocol 113, full inverse triangles). The full line is a guide 
for the eye. b) Normalized dielectric strength of the LF 
mode AeLF/(cDNA ■ Llf) ionic strength of the added salt 
(7s) for three representative DNA concentrations: cdna = 
0.05 mg/mL (diamonds, protocol III), cdna = 0.5 mg/mL 
(full triangles, protocol III; open triangles, protocol 112. 1; 
open inverse triangles, protocol 112. 2) and cdna = 0.8 mg/mL 
(circles, protocol 112. 1). 



IV. DISCUSSION 



chain. 



As a final remark, we point out that our results show- 
ing that the HF relaxation can be attributed to the 
semidilute mesh size ^, in the polyolcctrolyte semidilute 
solution in which single chain persistence length associ- 
ated with the LF relaxation is always larger than ^, con- 
firm the prediction of Odijk [19] that the same scaling 
law ^ oc c~°'^ should also be valid for semiflexible DNA 
polymers. 



First, let us summarize the results of dielectric spec- 
troscopy measurements. In the linear ac field regime, 
two broad relaxation modes are observed corresponding 
to three different time and length scales. The HF mode 
is centered in the frequency range between 0.1 MHz and 
15 MHz, depending solely on the DNA concentration as 
long as the DNA concentration remains larger than the 
added salt concentration. In this regime, the characteris- 
tic length scale is identified with the mesh size, and varies 
as ^ oc Cq^'^. Our data also seem to indicate that once 
the added salt becomes larger than the concentration of 
DNA intrinsic counterions, the high frequency charac- 
teristic length, rather than scaling as the semidilute so- 
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lution correlation length, levels off at a value close to 
the corresponding Debye length. More systematic exper- 
iments are needed to asses the possible added salt depen- 
dence in this regime of salt concentrations. In the limit of 
low DNA concentrations and low added salt, the semidi- 
lute solution correlation length smoothly crosses over to 
a less rapid scaling as oc c~^-^'^ probably reflecting the 
appearance of locally fluctuating regions with exposed 
hydrophobic cores. 

The LF mode is centered in the frequency range be- 
tween 0.5 kHz< z/LF <70 kHz. In DNA solutions with 
added salt, the characteristic length scale of this mode 
corresponds to the persistence length, which varies ex- 
perimentally . In the limit of low added salt, 
the characteristic length scale smoothly merges with the 
average size of the Gaussian chain composed of corre- 
lation blobs, which varies with DNA concentration as 

The dielectric data also seem to indicate that the free 
DNA counterions are mostly responsible for the high fre- 
quency relaxational mode, whereas the low frequency 
mode appears to be more complicated and the decou- 
pling of MO condensed and free DNA counterions seems 
to be difficult with any degree of confidence. 

Investigation of dielectric properties of DNA goes back 
to early '60s and since then a reasonable number of pa- 
pers have been published [24, 25, 30, 35, 43, 44, 45, 
46, 47]. In these studies, two dispersion modes associ- 
ated with the counterion fluctuations were found, one 
at very low frequencies, which depended on molecular 
weight (i.e. degree of polymerization, N) and another 
one in the intermediate frequency region, which showed 
no N dependence, but did show a pronounced concentra- 
tion dependence. For the sake of completeness we also 
mention the third relaxation in the GHz frequency range, 
which is due to solvent (water) relaxation and therefore 
not directly a concern of this paper. 

The low frequency mode was consistently associated 
with the DNA contour length, while the intermediate fre- 
quency one was often associated with the DNA statistical 
segment length. These interpretations were based on the 
theoretical models of the Mandel group [35, 36, 48] devel- 
oped for the case of a single polyelectrolyte (limit of very 
dilute solutions). Another interpretation of the inter- 
mediate frequency relaxation, proposed and verifled un- 
til now only for synthetic monodisperse polyelectrolytes, 
was that it was due to the counterion fluctuation along 
the correlation length [25, 41, 42]. It is worth noting the 
suggestion of Odijk [19] that the mesh size gives a phys- 
ical meaning to the statistically independent chain seg- 
ment length, which also scales as c^""'' [35]. The segment 
length was postulated to be due to the potential barriers 
along the chain, whose origin might be in the correlation 
length that measures the mean distance between contact 
points of the overlaping chains [20] . 

Our work reveals the existence of two, rather than 
one, distinct relaxation modes in the s.c. intermediate 
frequency range due to either a correlated response of 



counterions in the mesh of DNA chains in the solution, 
what we call the HF mode, or due to a response of the 
counterions along a single DNA chain, what we refer to 
as the LF mode. 

No previous experimental work was able to distinguish 
these two concentration dependent dispersions. The rea- 
son probably lies in the fact that the data analysis per- 
formed in the reported DNA dielectric spectroscopy stud- 
ies was not powerful enough to reveal and characterize 
two modes so close in frequency, where in addition one 
of them is small in amplitude. Indeed, work by Lee and 
Bone [47] hinted at two overlapping dispersions, but the 
authors were not able to characterize the smaller one 
properly. Another reason lies in the fact that none of 
these investigations covered so wide a range of DNA con- 
centrations and added salt ionic strengths, as we did in 
this work. 



A. Conformation of DNA in low salt solutions 

In this work we have paid special attention to the is- 
sue of the stability of ds-DNA helix, i.e. to the denatu- 
ration phase diagram [10]. The issue of DNA conforma- 
tion in pure water solutions is of paramount importance 
for proper understanding of our experiments. The ques- 
tion here is whether DNA at very low salt conditions 
is in the double stranded or single stranded form. Di- 
electric spectroscopy results strongly indicate that the 
double stranded form of DNA is stable in all pure water 
[27] solutions studied. Of course, a precise and defini- 
tive information on the polyelectrolyte intrachain confor- 
mation of DNA solutions would demand the small-angle 
neutron scattering and/or X-ray scattering experiments 
performed at the same conditions. 

First evidence for the stability of the ds-DNA form 
comes from the fact that the results of the dielectric 
spectroscopy measurements obtained in DNA solutions 
with added salt prepared from water (protocol III) co- 
incide with the results obtained on DNA solutions pre- 
pared from 10 mM NaCl (protocols 112, 113) (see Fig. 2b), 
3b); 4b) and 5b)). Second, we mention that Mandel [35] 
already reported that dielectric behavior in pure water 
DNA solutions was not found to differ markedly from 
that at low salt concentrations. Also the measured os- 
motic coefficient in nominally pure water conditions [49] 
confirms the assumption of an intact double stranded 
DNA form. Since the Manning charge density param- 
eters for ss- and ds-DNA are so different this difference 
should be apparent also in the measured osmotic coeffi- 
cient. None is detected however. 

In an attempt to clarify more this issue, we have mea- 
sured dielectric properties of pure water DNA solutions 
for the DNA concentration range between cdna = 0.5 
mg/mL and Cdna = 0.01 mg/mL, prepared according to 
protocol I, before and after the controlled denaturation 
protocol. Denaturation was accomplished by the heating 
of solutions for 20 min at temperature of 97°C, followed 
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by quenching to 4°C. Dielectric measurements were sub- 
sequently made at 25°C. The observed results were simi- 
lar for all studied DNA solutions. The dielectric strength 
and the relaxation time of the LF mode decreased sub- 
stantially after the heating, while the change observed for 
the HF mode was much smaller (Fig. 6). For cdna = 0.01 
mg/mL, the LF mode was not observed at all after the 
denaturation, implying the dielectric strength Ae < 1 
(see Fig. 6b)). Notice that Ae of the LF mode for this 
DNA concentration is about 10. 

We note that Llf, which measures the average size of 
the DNA chain, decreased after the heating, but showed 
the same power law behavior Llf oc Cqna^^" 
the untreated DNA solution (Fig. 7b)). The observed 
change indicates that denatured ss-DNA, which is in the 
form of coil, is shorter and has smaller average size of the 
chain. Furthermore, -Lhf which measures the correlation 
length of the DNA solution, decreased after the heating 
and showed the power law behavior Lrf occ^^'^'^ in the 
whole DNA concentration range (Fig. 7a)). A smaller 
value of LiiF is expected for the denatured DNA solu- 
tion, since such a solution should contain twice the num- 
ber chains. The power law behavior with the exponent 
-0.33, observed now also for the larger DNA concentra- 
tions, indicates that the hydrophobic core of DNA is fully 
exposed once the heating protocol was applied. 

All these results confirm that although untreated 
DNA at low salt and semidilute conditions might 
show locally exposed hydrophobic cores in a dynamic 
sense, a real/complete unzipping and separation of the 
strands might be accomplished only after the denatu- 
ration/heating protocol is applied. These observations 
therefore suggest that DNA in pure water solutions is 
not denatured into two spatially well separated single 
strands, but is rather in the double stranded form, lo- 
cally interspersed with exposed hydrophobic cores in the 
limit of low DNA concentrations. 

Furthermore, even after denaturation the two DNA 
strands appear to remain in relatively close proximity 
rather than becoming completely dissociated, an obser- 
vation well substantiated also by the correlation length 
measurements by SANS at semidilute DNA conditions 
[50]. In these measurements the correlation length mea- 
sures the characteristic distance between the hydrogen- 
containing (sugar-amine base) groups. Hammouda and 
Worcester [50] have recently determined that on melting 
of DNA in DNA/d-ethylene-glycol mixtures the correla- 
tion length increases from about 8 A to about 12-15A, 
which implies that even after melting the two strands of 
the ds-DNA remain in very close proximity. This appears 
to be due to the presence of other chains in the semidi- 
lute solution that spatially constrain the separate strands 
even after melting and prevent complete dissocation [51]. 

The question furthermore arises if the ronaturation 
process is fast enough to occur partially during dielectric 
measurements which are performed at 25°C. In order to 
check this, we have repeatedly measured the response of 
pure water DNA solution, after the heating protocol to 



100 



=Vj 10 r 



1 



a) 



'"I I I I 

pure water DNA solution 

25 °C 



unti-eated 




denatured at: \ 



I I /l A I 



101 102 10^ 10^ 105 10^ 10' 108 
V(Hz) 



100 



10 r 



■^0 



: 1 1 

;b) 


1 1 1 1 

pure water DNA solution 




25 °C 


■ untreated 








denatured at: 

■ 97°C 

1 1 . . 





0.1 

101 102 103 104 10^ 106 10' 10» 
V(Hz) 

FIG. 6: Double logarithmic plot of the frequency depen- 
dence of the imaginary part of the dielectric function (e" ) at 
T=25°C of a pure water DNA solution (protocol I) with DNA 
concentration cdna = 0.1 mg/mL (panel a)) and cdna = 0.01 
mg/mL (panel b)) before (denoted as untreated) and after the 
heating to 97°C (denoted as denaturated, see Text). The full 
lines are fits to the sum of the two HN forms; the dashed lines 
represent a single HN form. 



97 "C was applied. The time span was 100 min, the first 

measurement was taken 3 minutes after the sample was 
heated from 4°C to 25 °C. The observed response after 3 
min and after 100 min was the same (inside the error bar 
of 1.5%) indicating very long time constant characteriz- 
ing the renaturation process. This result confirmed that 
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FIG. 7: a) The characteristic length Lhf of the HF mode 
before (open squares) and after the heating to 97° C (fuU 
squares). The full line is a fit to the power law Lhf ocCq^'^^ 
and ocCq^'^ for cdna smaller and larger than Cco ~ 0.6 
mg/mL, respectively. The dashed line is a fit to the power 
law Lhf ocCqJ^I,'^'* in the whole DNA concentration range, 
b) The characteristic length Llf of the LF mode before 
(open triangles) and after the heating to 97°C (full trian- 
gles). The full line and dashed lines are fits to the power law 

r „-0.29±0.04 
i^LF OC CjjNA 



ds-DNA was found to be stable at 25°C dissolved in nom- 
inally pure water (no added salt, 0) for the concen- 
tration of intrinsic counterions larger than 0.2 mM. This 
again is consistent with our experiments on dielectric re- 
laxation of DNA solutions. 

Another meaningful question would be, what is the 
smallest c,„ which can still keep DNA in the double 
stranded form. The UV absorbance results indicated 
that no added salt is needed to stabilize ds-DNA in the 
case when the intrinsic counterion concentration, c,„, 
is larger than 0.2 mM, while dielectric spectroscopy re- 
sults suggest it can be one order of magnitude smaller, 
Cin > 0.03 mM. Notice that in the latter CclSG C-iri IS still 
larger than the estimated ion concentration of pure wa- 
ter (2/s = 0.02 mM), i.e. Cin > 21 s- The osmotic pres- 
sure data indicate [52] that intrinsic counterions them- 
selves can stabilize ds-DNA for Cm larger than 30 mM 
and 150 mM for two different added salt concentrations 
27s = 4 mM and 20 mM, respectively. All these results 
suggest that the concentrations of intrinsic counterions 
and added salt themselves are less important, rather it 
is their ratio that defines which limit prevails. Thus ab- 
sorbance and osmotic pressure data both suggest that 
for Cj„/2/s > 10 intrinsic counterions prevail in stabiliz- 
ing ds-DNA helix, while dielectric results suggest that 
this ratio might be shifted to Ci„/2/s > 1. Neverthe- 
less, the unusual scaling exponent -0.33 of the semidilute 
correlation length found for the DNA pure water solu- 
tions with DNA concentrations smaller than Cco ~ 0.6 
mg/mL seem to suggest the existence of local conforma- 
tional fluctuations which partially expose the hydropho- 
bic core of DNA. In order to verify this proposal a further 
comparative study, including UV spectrophotometry and 
dielectric spectroscopy, is planned to study in depth the 
conformation of DNA in pure water solutions. We are 
however convinced that in all our investigations DNA 
did not exist as two separated single coils, rather DNA 
was effectively always in its double stranded form. 



the observed dielectric properties are the ones of ss-DNA. 

Osmotic pressure resiflts of Raspaud et al. [52] on short 
nucleosomal fragments of DNA also indicate that at high 
enough concentration, intrinsic DNA counterions prevent 
destabilization of ds-DNA helix in pure water solutions. 
These experiments show that the role of added salt {Is) in 
the ds-DNA stabilization becomes negligible in the high 
DNA concentration range i.e. in the limit Cdna ^ Is- 
Osmotic pressure of intrinsic DNA counterions can thus 
be high enough to prevent denaturation of ds-DNA helix 
in pure water solutions. 

UV spectrophotometry experiments performed previ- 
ously by Record [53] on T4 and T7 phage DNA have 
also shown that ds-DNA denaturation depends not only 
on added salt concentration 7^, but also on the con- 
centration of intrinsic DNA counterions, Ci„ [mM] = 
Cdna [mg/mL] x 3 /xmol/mg (see Section 11). Moreover, 



B. Ionic screening: added salt versus intrinsic 
DNA counterions 

Next we address the issue of the respective roles in 
ionic screening of intrinsic DNA counterions and ions 
from the added salt. First, we examine conditions un- 
der which the OSF expression for the persistence length 
Lp = Lq + a ■ I^^ is valid. Dielectric data show that 
the influence of the added salt on the persistence length 
is important as long as the ionic strength 7^ is suffi- 
ciently larger than the concentration of intrinsic coun- 
terions. Plot in Fig. 8 reveals that this condition reads 
27s > O.Acin- In this limit where the OSF theory applies, 
the coefficient a should be equal 0.324 AM (assuming 
MO counterion condensation). However, experimentally, 
different values are found. Measurements of DNA elastic 
properties as a function of ionic strength also yielded the 
coefficient in the OSF expression, a = 0.8 Am [17], which 
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is different from the one expected from the MO counte- 
rion condensation theory. On the other hand, a magnetic 
birefringence study by Maret and Weill [54] suggested the 
values of a between 0.25 and 0.45, therefore indeed close 
to 0.324. However, the authors fitted their data to the 
(0.12ci„+/s)~^ instead to the OSF I^^ dependence. The 
first term represents the influence of free DNA counteri- 
ons. It is rather obvious that the OSF fit would not yield 
the MO value of the coefficient a. Why thus is there a dif- 
ference between the MO value and the measured value of 
the coefficient in the OSF dependence of the persistence 
length on the ionic strength of the added salt? 

In an attempt to reconcile these various values, we have 
rescaled the persistence length from dielectric measure- 
ments as 2.5 • Lup — > Llf in order to collapse the behav- 
ior from dielectric and elastic experiments onto a single 
curve (Fig. 9). Rescaling is justified since the expres- 
sion connecting Lp from dielectric properties with the 
measured mean relaxation time is only valid as a scaling 
relationship. Numerical factors are less straightforward 
and are essentially unknown. An exact expression to- 
gether with an appropriate numerical coefficient, known 
as the Einstein- Smohichowski formula, is only known in 
the dc limit, where L = \J2Dt. 

The fit of this new rescaled version of the Llf to 
the OSF expression with two free parameters now gave 
Lo = 530 A and a = 0.69 AM. The difference between 
this value of a and the one expected from the MO counte- 
rion condensation theory signals a different effective lin- 
ear charge density than the one stemming from the MO 
countcrion condensation theory. This difference in the 
measured and theoretically expected value of a coefficient 
of the OSF fitting form for the persistence length would 
be consistent also with other experiments that provide a 
value for the effective charge density [14]. 

Finally, let us examine the low salt limit (2/^ < OAcin) 
in which we expect that the intrinsic counterions be- 
come dominant. The characteristic length of the LF 
mode varies with DNA concentration as Xlf oc c^^^, 
or equivalently as c^°'^^. This value of the exponent 
would allow us to identify _Llf with the average size of 
the Gaussian chain composed of correlation blobs scaling 
as Cqna^' where electrostatic interactions are screened by 
other chains and counterions, and thus DNA acts as its 
own salt. 



V. SUMMARY 

All these results make it quite clear that for the dielec- 
tric properties studied in this contribution the role of free 
and condensed intrinsic DNA counterions can not always 
be separated. For the HF relaxation we are reasonably 
sure that it is mostly the free counterions that act as the 
relaxation entities in a type of hopping process between 
the correlated DNA chains in a semidilute polyelectrolyte 
solution. On the contrary, for the LF relaxation process 
it seems to be mostly the condensed counterions that re- 
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FIG. 8: Characteristic length of the LF mode (Llf) nor- 
malized with tiie value in pure water solutions (low salt limit 
Is ~ 0.01 mM) vs. added salt concentration normalized by 
the concentration of intrinsic counterions {2Is/cin)- Data 
are for three representative DNA concentrations: cdna = 
0.05 mg/mL (diamonds, protocol III), cdna = 0.5 mg/mL 
(full triangles, protocol III; open triangles, protocol 112. 1; 
open inverse triangles, protocol 112. 2) and cdna = 0.8 mg/mL 
(circles, protocol 112. 1). 
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FIG. 9: Characteristic length of the LF mode, rescaled as 
2.5-Llf — > ^LF {LuF data are from Fig. 4b)) is shown together 
with data from pulling of DNA (crosses) taken from Ref. [17]. 
The full line and dashed lines are fits to OSF theory with 
Lo = 530 A and a = 0.69 AM and with Lo = 500 A and a = 
0.324 Am corresponding to the theoretical OSF prediction 
with complete MO condensation. 



lax along or in close proximity to an individual chain in 

the polyelectrolyte solution, except at large enough salt 
concentrations where at least some of the free counterions 
seem to carry on this role. 

For the HF relaxation, the experimentally found spa- 
tial correlation providing the characteristic size of the 
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relaxation process is given universally by the polyelec- 
trolyte semidilute solution correlation length. The LF 
relaxation can be characterized by two different correla- 
tion scales. At low salt concentration this correlation 
length can be identified with the average size of the 
polyelcctrolyte chain, where electrostatic interactions are 
screened by other chains and counterions, thus where 
DNA acts as its own salt. At higher salt concentra- 
tions, the spatial correlations arc provided by the single 
chain orientational correlation length, i. e. the persistence 
length, that depends on the salt concentration via the 
OSF mechanism. Both the correlation lengths for the 
LF relaxation seem to be telling us, that it is the single 
chain that is responsible for the relaxation process, in 
which counterions either fluctuate orientationally on the 
length scale of the persistence length, or they fluctuate 
between two ends of the average size of the chain in the 
solution. 

In conclusion, our results demonstrate that there are 

three fundamental length scales that determine the di- 
electric response of a semidilute DNA solution: the av- 
erage size of the polyelectrolyte chain in the regime in 
which DNA acts as its own salt, the OSF salt-dependent 
persistence length of a single polyelectrolyte chain and 
the dGPD semidilute solution correlation length or the 
mesh size of the polyelectrolyte solution. While the free 
DNA counterions can be identified reasonably well as 
the relaxation entities of the HF relaxation mode, the 
LF relaxation mode does not allow for such a clearcut 



separation between MO condensed and free counterions 
as relaxation entities. Our data suggest that in fact 
both of them contribute to various extent in different 
salt and DNA concentration regimes. Finally, our results 
confirm that although double-stranded DNA at low salt 
concentrations shows locally exposed hydrophobic cores 
in a dynamic sense, unzipping of the strands is accom- 
plished only after the denaturation/heating protocol is 
applied. But even then this unzipping of the two DNA 
strands is probably at most local and complete separa- 
tion of the strands at semidilute solutions is never really 
accomplished. This issue will be addressed in our further 
studies. 



Acknowledgments 

ICP-AES was performed by M. Ujcvic at the Croa- 
tian National Institute of Public Health. We thank 
T. Ivek, D. Vurnek and R. Zaja for help in the data 
analysis and electrophoresis measurements. Discussions 
with D. Baigl, F. Livolant and E. Raspaud are greatly ac- 
knowledged. This work was supported by the Croatian 
Ministry of Science, Education and Sports. Rudi Pod- 
gornik would like to acknowledge the financial support 
of the Agency for Research and Development of Slovenia 
under grant P1-0055(C). 



[1] M. Daune, Molecular Biophysics (Oxford University 

Press, New York, 2003). 
[2] F. Oosawa, Poly electrolytes, (Marcel Dekker, New York, 

1971). 

[3] K. S. Schmitz, Macroions in Solution and Colloidal Sus- 
pension (VCH New York, 1993). 
[4] A. V. Dobrynin and M. Rubinstein, Frog. Polym. Sci. 30, 
1049 (2005); A. V. Dobrynin, R. H. Colby and M. Ru- 
binstein, Macromolecules 28, 1859 (1995). 
[5] R. R. Netz and D. Andelman, Pliys. Rep. 380, 1 (2003). 
[6] H. Boroudjcrdi, Y. -W. Kim, A. Naji, R. R. Netz, 
X. Schlagberger and A. Serr, Phys. Rep. 416, 129 (2005). 
[7] P. L. Hansen, D. Svensek, V. A. Parsegian and R. Pod- 

gornik, Phys. Rev. E 60, 1956 (1999). 
[8] A. Y. Grosberg, T. T. Nguyen and B. I. Shklovskii, 

Rev. Mod. Phys. 74, 329 (2002). 
[9] P. L. Hansen, R. Podgornik and V. A. Parsegian, 
Phys. Rev. E 64, 021907 (2001). 
[10] V. A. Bloomfield, D. M. Crothers and I. Tiuocco, Jr., Nu- 
cleic Acids (University Science Books, Sausalito, 2000). 
[11] R. Das, T. T. Mills, L. W. Kwok, G. S. Maskel, I. S. Mil- 
Ictt, S. Doniach, K. D. Finkelstoin, D. Herschlag and 
L. Pollack, Phys. Rev. Lett. 90, 188103 (2003). 
[12] H. H. Strcy, V. A. Parsegian and R. Podgornik, 

Phys. Rev. E 59, 999 (1999). 
[13] R. Podgornik, D. C. Rau and V. A. Parsegian, Macro- 
molecules 22, 1780 (1989). 
[14] R. Podgornik, D. C. Rau and V. A. Parsegian, Bio- 



phys. J. 66, 962 (1994). 
[15] T. Odijk, J. Polym. Sci.: Polym. Phys. 15, 477 (1977). 
[16] J. Skolnick and M. Fixman, Macromolecules 10, 944 

(1977). 

[17] C. G. Bauinaim, S. B. Smith, V. A. Bloomfield and 

C. Bustamante, Proc. Natl. Acad. Sci. USA 94, 6185 
(1997). 

[18] M. UUner, J. Phys. Chem. B 107, 8097 (2003). 

[19] T. Odijk, Macromolecules 12, 688 (1979) 

[20] P. G. do Gcnrics, P. Pincus, R. M. Velasco and 
F. Brochard, J. Phys. (Paris) 37, 1461 (1976). 

[21] P. Pfeuty, J. Phys. (Paris) 39, C2-149 (1978). 

[22] C. F. Schmidt, M. Barmann, G. Isenberg and E. Sack- 
mann, Macromolecules 22, 3638 (1989). 

[23] R. Pethig, Dielectric and Electronic Properties of Biolog- 
ical Materials (J. Wiley, New York, 1979). 

[24] N. Naudi, K. Bhattacharyya and B. Bagchi, 
Chem. Rev. 100, 2013 (2000). 

[25] F. Bordi, C. Camctti and R. H. Colby, J. Phys.: Con- 
dens. Matter 16, R1423 (2004). 

[26] S. Tomic, T. Vuletic, S. Dolanski Babic, S. Krca, 

D. Ivankovic, L. Griparic and R. Podgornik, 
Phys. Rev. Lett. 97, 098303 (2006). 

[27] Throughout this article by pure water we mean 
ultrapure MilliQ water with declared con- 
ductivity a-=0.056 fjS/cm; http://www. milli- 
pore. com/markets/bioscienco. nsf/home. 

[28] R. F. Steiner and R. F. Beers, Polynucleotides (American 



15 



Elsevier, New York, 1961). 
[29] W. Beers, A. Ccrami, and E. Reich, 

Proc. Natl. Acad. Sci. USA 58, 1624 (1967). 
[30] B. Sail, R. K. Mohr, C. J. Montrose and T. A. Litovitz, 

Biopolymers 31, 1171 (1991). 
[31] S. Havriliak and S. Negami, J. Polym. Sci. C 14, 99 

(1966). 

[32] M. Pintoric, T. Vulctic, S. Tomic and J. U. von Schiitz, 

Eur. Phys. J. B 22, 335 (2001). 
[33] R. Hinrichs and J. A. H. da Jornada, Rev. Sci. In- 

strum. 68, 193 (1997). 
[34] W. E. Sonnen, G. E. Wesemberg and W. E. Vaughan, 

Biophys. Chem. 32, 283 (1988). 
[35] M. Mandel, Ann. NY Acad. Sci. 303, 74 (1977). 
[36] M. Mandel and T. Odijk, Ann. Rev. Phys. Chem. 35, 75 

(1984). 

[37] S. S. Dukhin and V. N. Shilov, Adv. Coll. Interface 
Sci. 13, 153 (1980). 

[38] R. W. O'Brian, J. Coll. Interface Sci. 113, 81 (1986). 

[39] In the literature intrinsic DNA counterions arc classified 
into two groups bearing different names. The first is com- 
monly named as condensed or bound and the second as 
diffuse or free. Throughout this article we will refer to 
these two groups as condensed and free counterions. 

[40] T. E. AngeUni, R. Golestanian, R. H. Coridan, J. C. But- 
ler, A. Beraud, M. Krisch, H. Sinn, K. S. Schweizer and 
G. C. L. Wong, Proc. Natl. Acad. Sci. USA 103, 7962 
(2006). 

[41] K. Ito, A. Yagi, N. Ookubo and R. Hayakawa, Macro- 



molecules 23, 857 (1990). 
[42] F. Bordi, C. Carnetti, T. Gili and R. H. Colby, Langrrmir 

18, 6404 (2002). 
[43] M. Sakamoto, H. Kanda, R. Hayakawa and Y. Wada, 

BiopolymerslS, 879 (1976). 
[44] R. J. Molinari, R. H. Cole and J. H. Gibbs, Biopolymers 

20, 977 (1981). 

[45] S. Takashima, C. Gabriel, R. J. Sheppard and 
E. H. Grant, Biophysical J. 46, 29 (1984). 

[46] S. Bone and C. A. Small, Biochim. Biophys. Acta 1260, 
85 (1995). 

[47] R. Lee and S. Bone, Biochim. Biophys. Acta 1397, 316 
(1998). 

[48] F. Van der Touw and M. Mandel, Biophys. Chem. 2, 218 

(1974). 

[49] H. E. Auer and Z. Alexandrowicz, Biopolymers 8, 1 
(1969). 

[50] B. Hammouda and D. Worcester, Biophys. J. 91, 2237 
(2006). 

[51] B. Hammouda, personal communication (2006). 

[52] E. Raspaud, M. da Conceiao and F. Livolant, 

Phys. Rev. Lett. 84, 2533 (2000). 
[53] M. T. Record, Jr., Biopolymers 14, 2137 (1975). 
[54] G. Maret and G. Weill, Biopolymers 22, 2727 (1983). 
[55] M. N. Spiteri, F. Boue, A. Lapp and J. P. Cotton, 

Phys. Rev. Lett. 77, 5218 (1996). 
[56] M. Muthukumar, J. Chem. Phys. 105, 5183 (1996). 



